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Abstract 

The degenerate Anderson model with a nonlinear electron dispersion and an 
energy dependent hybridization is proven to exhibit hidden integrability and is 
diagonalized by the Bethe ansatz. If the impurity /-level energy lies below the 
conduction band edge, the spectrum of the system is shown to contain pairs 
of confined "gap excitations" and their heavy bound complexes. In addition, 
multielectron localized states due to a pinning of gap bound complexes to the 
impurity are predicted. 
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It is well known that the basic theoretical models describing magnetic impurities in non- 
magnetic metals [JJ, such as the s-d (Kondo) model, the Anderson model, etc., are integrable 
under two additional conditions: (i) an electron-impurity coupling is assumed to be energy 
independent, and (ii) a band electron dispersion e(k) can be linearized around the Fermi 
level, e(k) ~ vp(k — kp), where kp and Vp are the Fermi momentum and velocity, respec- 
tively. Only under these conditions, both an electron-impurity scattering and an effective 
electron-electron coupling are described in terms of discontinuous jumps in the Bethe ansatz 
wave functions. Therefore a linear particle spectrum and a pointlike particle-impurity cou- 
pling are considered now as the necessary mathematical conditions for integrability of the 
"impurity" models. In what follows we will use the phrase "linear dispersion approximation" 
(LDA) to refer to a model modified in accordance with these conditions. The LDA is clear to 
eliminate both the effects of the band structure and possible excitations of the system with 
energies lying inside the interband gap from an exact analysis of the behaviour of magnetic 
alloys. Therefore an extension of the Bethe ansatz technique to models accounting for a 
more realistic band structure and an energy dependence of particle-impurity coupling is of 
great physical interest. 

In the present paper, we study the n-fold degenerate Anderson model with a nonlinear 
band electron dispersion and an energy dependent electron-impurity coupling, which is the 
basic theoretical model for mixed- valent, Kondo and heavy-fermion phenomena JE]. The 
model Hamiltonian is written in the following one-dimensional form: 

} + e/^X aQ . (la) 

The Fermi operator c^ a {k) creates a conduction band electron with the total angular momen- 
tum projection a and the momentum modulus k. A rare-earth impurity is described by the 
Hubbard operators X ab with algebra X ab X cd = 6 cb X ad , where the index a = 0, a enumerates 
both the nonmagnetic state (a = 0) and the degenerate magnetic states (a = 1, . . . ,n) of 
the impurity. The third term is the /-level energy of an electron localized on the impurity, 
while the second one describes the hybridization of band and /-level electron states with 
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— \e{k)ci{k)c a {k) + v{k) cl(k)X 0a + X a0 c a (k) 



the matrix element v(k). The Coulomb repulsion on the impurity orbital is assumed to be 
very large, so that the multiple occupancy of the impurity /-level is excluded. Within LDA 
the Anderson model has been diagonalized by Wiegmann |2j (nondegenerate version) and 
by Schlottmann || (degenerate version). 

It has recently been discerned || that LDA is not necessary for an exact diagonalization 
of the basic impurity models of quantum optics, describing a system of Bose particles with 
a nonlinear dispersion coupled to two-level atoms. Making use of some mathematical analo- 
gies between "magnetic" and "optical" models, we diagonalize exactly the Hamiltonian (la) 
and derive a hierarchy of the Bethe ansatz equations. The information about the nonlin- 
ear electron dispersion and the energy dependence of the hybridization is contained only in 
the nonlinear energy dependence of "rapidities" of "charge" excitations of the system. If 
the impurity /-level energy £f lies below the conduction band edge, the spectrum of the 
system exhibits a novel extremely rich class of charge excitations with energies lying inside 
the interband gap around 6f. These excitations have been excluded obviously from previous 
studies by LDA. The gap spectrum is shown to contain propagating pairs of confined "gap 
excitations", which do not exist separately from each other, and their heavy bound com- 
plexes. In addition to propagating complexes, we demonstrate the existence of multielectron 
gap states localized in the vicinity of the impurity due to a pinning of a gap bound complex 
to the impurity. 

To diagonalize the model Hamiltonian, it is convenient first to rewrite Eq. (la) in 
terms of electron operators on the "energy scale", c a (e) = (de(k) / ' dk)" 1 ^ ] c a (k) , with algebra 
{c a (e), cl,(e')} = 2ir6 aa >S(e - e') 7 

} + e/ ^X aa , (lb) 

a 

where V(e) = (de(k) / dk)~ 1 ^ 2 v(k) . Note that the energy dependence of the effective coupling 
V(e) is determined now by both the energy dependence of the hybridization matrix element 
and the nonlinearity of the electron dispersion. 



Ucl(e)c a (e)V(e) cl(e)X 0a + X a0 c a (e) 
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Now we look for one-particle eigenstates of the model in the form 



l*i> = E^ 



POO rip 

5*«o + 1 ^V{e)4>{e)ci{e) 



10), 



where A a is arbitrary, and the vacuum state contains no band electrons and the impurity is 
in the nonmagnetic state, c a (e)|0) = X aa \0) = 0. The Schrodinger equation then reads 

{-id T -u)<l>(T\u) + g(u)5(T)=0, (2a) 

/oo 
drV{T)4>{r\u) = 0, (2b) 
-oo 

where 0(r) is the Fourier image of the auxiliary wave function <f)(e), while the effective 
particle- impurity coupling T(r) contains the information about the electron dispersion, 

de % t — i / \ f 00 de. 



T 



^-4>{e)e l£ \ r(r) = / ^V\e)e 
oo Ztt Jo /vt 



For uj > 0, one gets g(a>) = + z/2] 1 and 

MM, A ~ (V 2 ) s g n ( r ) 

where the "rapidity" 

h(u) = (u + X'{u)-e f )/V 2 (u), 

while S'(o;) and £"(ti>) = V r2 (cj)/2 are respectively the real and imaginary parts of the self- 
energy £(to>) = J °°(de/27!-)V 2 (e)/ (e — uj — iO). In what follows we neglect the uj dependence 
of E'(o;), S'(o;) — > £'(e/), and regard £'( e /) as a small correction to the /-level energy due 
to the hybridization with the band electron states, ej = tf — Y!(oj). If 6f lies below the band 
edge, 6f < 0, Eqs. (2) admit also the discrete mode with the eigenenergy ujj, = e/, which 
obviously corresponds to the electron-impurity bound state. In the case of a Bose system 
with gap dispersion, a particle-impurity bound state has earlier been predicted by John and 
Wang ||. To keep notations simple, hereafter the bar in e/ is omitted. 

The auxiliary wave function (3) is discontinuous, but the electron wave function 
ip a (x\u>) = (c a (x)\^x) = i])(x\uj)A a (where the electron operator c a (x) is defined as 
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c a (x) = f£°(dk/2Ti)c a (k) exp (ikx) [JJ) is continuous and results from the integral "dress- 
ing" of the auxiliary function, 



ip(x\u) = / dru(x; t)(/)(t\u), 

J —oo 

with the dressing function 

_.n.\ ( d s(k) \ i(kx~e{k)r) 



u(x;t) = J 7^ v ( k ) 



/o 2tt \ dk J 

In LDA, the dressing function is nothing but the delta-function, u(x;t) ~ S(x — r), and 
hence the auxiliary wave function <f>(r\u) is the same as the electron wave function ip(x\u). 

The idea of auxiliary functions (or "auxiliary particles") can easily be extended to the 
multielectron case. For instance, the two-electron wave functions in the energy space are 
represented as ^ ai a 2 { £ ii e 2) — V( £ i)V( £ 2)$a 1 a 2 (. £ i, e 2) an d 4 lQ2 (e) = V(e)G aia2 (e), where 
the latter describes the state, in which one of the electrons is localized on the impurity. In the 
auxiliary r-space, the Schrodinger equations for the Fourier images of the auxiliary functions 
$ Q , 1Q , 2 (r 1 , r 2 ) and G aia2 (r) are then similar to those within LDA but with the nonlocal 
particle-impurity coupling T(t). It is remarkable that, despite the nonlocal coupling, the 
two-particle scattering matrix is still found ||] in the well-known form: 

g = K^i) - h(u 2 )+iP (A) 
h{u x ) - h(uj 2 )+i 

where P = 5 aia i8 a20l i is the permutation operator. The ^-matrix satisfies obviously the 
Yang-Baxter equations. Hence the multiparticle scattering of the auxiliary particles is fac- 
torized into two-particle ones and the auxiliary multiparticle wave functions have the ordi- 
nary Bethe ansatz structure. But due to the nonlinear electron dispersion and the nonlocal 
electron-impurity coupling, the multielectron wave functions 

r-oo N 

ty ai ...a N (xi, . . . ,X N ) = / ® ai ...a 2 (n, ■ ■ ■ ,T N ) ]J u(xf, T^dTj (5) 

J-oc j=1 

are continuous functions of the coordinates Xj. The factorization of multielectron scattering 
and the Bethe ansatz construction for multielectron wave functions are thus hidden and 
manifested only in the limit of large interelectron separations. 



Imposing the periodic boundary conditions on the A-electron wave function (5) on the 
interval of size L, we arrive at the following Bethe ansatz equations (BAE): 

Mi 

exp {ikjL)e x {\f ) = n ei(Af - A«) (6a) 

ai=l 

Mi+„ Mi 

n n e i(AS-A^) = -ne 2 (A2-A?) (6b) 

iv=±l a; + „=l 6=1 

where £/ = Sj^j is the eigenenergy, kj = k(ujj) is the momentum of a charge excitation 
with the energy Uj, and e n (x) = (x — in/2)/(x + in/2). If is the number of particles 
with the color a, the numbers Mi are defined by Mi = J2a=} ™ai being the total number 
of electrons, M = AT. It is clear that only charge excitations with rapidities = ^(cUj) 
contain the information about the electron dispersion and the energy dependence of the 
hybridization, while BAE for the "color" rapidities {A®}, I = 1, . . . , n — 1 coincide with the 
corresponding equations in LDA. Therefore, we can use previous results || to classify the 
solutions of Eqs. (6) as L — > oo and we may focus our studies on charge excitations only. 

The charge rapidities can be grouped into a "string", which describe a bound complex 
of electrons with different colors, 

h j = A^+ l -{m + l-2j) ) (7) 

where j = 1, ... ,1 + 1 enumerates now the string rapidities, and an arbitrary real A^ denotes 
a rapidity of corresponding color excitation. Since only electrons of different colors interact, 
one can build strings containing up to m = I + 1 < n charge excitations. Eq. (7) is a 
solution of BAE, if and only if the imaginary parts of the rapidity hj and the corresponding 
momentum kj have the same sign, 

sgn (Im hj) = sgn (Im kj) . (8) 

It is easy to understand that solutions of Eq. (8) determine vicinities on the u axis with an 
attractive interaction between charge excitations. In the other vicinities, where NC is not 
satisfied, the charge excitation coupling is repulsive, and hence they comprise no strings, 
except one-particle ones. Thus, though BAE are valid at arbitrary e(k) and v(k), their 
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solutions are controlled by the properties of analytical continuations of the functions k{ui) 
and h(u>) in the complex u;-plane. In LDA, where kj = uij and hj = (cUj — ej)/v 2 (fc F ), 
the necessary condition (NC) (8) is obviously satisfied at all frequencies. For the sake 
of simplicity, we confine ourselves to the case of a weak momentum dependence of the 
hybridization matrix element, assuming that v(k) can be replaced by some constant value 
y/2v. Moreover, we also assume that the electron dispersion has the simple form e(k) = k 2 , 
where we set the effective electron mass equals to 1/2. Then one obtains 

k{uS) = y/cu, h{uj) = v~ 2 {uj — £f)k{u)). (9) 

The analytical continuations of the functions k(u) and h{ui) depend essentially on the 
position of the real part of the frequency with respect to the band edge. If the real parts 
of all the frequencies uij corresponding to the /i-string (7) are positive, BAE with Eq. (9) 
lead only to some corrections to the system spectrum found within LDA. For the states of 
physical interest near the Fermi energy, these corrections are negligibly small. 

But the situation is changed drastically, if the real parts of all the frequencies ujj are 
assumed to lie below the band edge, Uj < 0. Let uo = £ + irj, where £ < 0, and k(£ ± 
iO) = ±i/t(£), where «;(£) = To avoid tedious algebraic expressions, we restrict 

further analysis to the case rj <C £. Then k(u) ~ sgn(r/)[— t]k'{^) + and h{u) ~ 

sgn(r])[— r]r'(£) + ir(£)], where r(£) = f~ 2 (£ — e/)/c(£). Since /c(£) is positive, NC leads now 
to the condition r(£) > 0, which is met only if e/ is negative and £ G (e/,0). Hereafter, 
we study only gap states lying around e/, because they are well separated from the band 
edge. The remaining analysis is then simplified by linearizing the function r(£) at the point 

f = e /> r (0 = a (£ ~ e /)> where a = v ~ 2K ( e f)- 

Let us start with the simplest case of a two-particle string. Its complex conjugated 

rapidities, hi = + i/2 and h 2 — h\ — A^ — i/2 are mapped to the corresponding 

pair of frequencies, uj = £ + irj and u* = £ — it], where rj is chosen to be positive, r\ > 0. 

The pair of corresponding momenta are given by k = q + and k* = q — where 

q = —t]k'{^) = t)\k'(£)\, because «;'(£) is negative. The real and imaginary parts of frequencies 
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are easily found as £ = e/ + l/2a and 77 = —A^/a. Since 77 > 0, only strings with 
negative A* 1 ) are mapped to the gap states. For the gap states the real part of rapidities 
A* 1 ) determines the imaginary part of frequencies rj, therefore strings with different A^ are 
mapped to pairs with the same real part £ but different imaginary parts rjp = A^/a. 

The expressions obtained describe a confined state of two gap excitations of the system, 
because, unlike excitations bound into ordinary band two-particle complexes, gap excitations 
do not exist separately from each other. These results are easily extended to the case of 
strings containing an even number of particles, m = 2/i. All pairs of complex conjugated 
rapidities, hj and h*, j = 1, ...,//, are mapped to the corresponding pairs of complex 
conjugated frequencies, uj = ^ + ir] and uo* = ^ — ir], and momenta kj = qj + 
and k* = qj - i^f). Here = e/ + (// + 1/2 - j)/a, 77 = A® /a and qj = t]\k'(^ 0) )\. 
The expressions obtained describe obviously the bound complex of \x pairs of confined gap 
excitations with the eigenenergy per particle 

4 0) = 2Eef = e/ + /i/2a. (10) 
i=i 

Note that the frequencies and momenta of the gap complex do not have the string structure. 
The upper index in the above expressions, (0), indicates that they are derived with linearizing 
of the function r(£) at the point £ = e/. In this approximation, the energy of a complex is 
independent of its momentum per particle, 

Q^^itq^ivi/riitwtf)], (11) 

and hence the states obtained are infinitely degenerate. But the degeneracy is lifted by 
the first correction to Eq. (10). To find this correction, one has to keep the next term, 
— i(?7 2 /2)r"(£), of the Taylor series for the function h(uj) and the next term in the expansion 
of the function r(£) at the point £ = ef. r(£) ~ a(£ — ej) + 6(£ — e^) 2 , where b = v~ 2 n'(ef) < 0. 
We then find the first order correction to the energy of a motionless complex and determine 
the kinetic energy contribution to the total energy of a complex in the effective mass ap- 
proximation, 
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E/i = — 8^ + Q 2 /2m* fl . (12) 

Here 5^ = j^wi^ 2 - 1) and m* = ^ Ej (K(£ (0) )l) 2 are the band half-width and the 
effective mass. For small /x, the bands of propagating complex are very narrow and complexes 
are very heavy and even motionless at Q — 0. But the bandwidth grows as /i 2 , so that large 
fj, complexes are quite mobile for large Q. At arbitrary Q, one needs to study the exact 
equations for the parameters of a complex, Re h{£j, rjj) = A"> and Im h(^j, rjj) = (ji—j+1/2). 

In the case of an /Vstring containing an odd number of particles, m = 2/i + 1, the extra 
real rapidity is mapped only to the gap state with £ = u d = ej, which corresponds to the 
electron-impurity bound state. The remaining complex conjugated pairs of rapidities of an 
odd string with AW — > + are then mapped to a deformed motionless gap complex with 
£j = £f + ((/, — j + 1). Since one of the particles of this complex is bound to the impurity, 
a complex as a whole is pinned to the impurity. A pinned complex describes obviously the 
state of the system, in which electrons with different colors are localized in the vicinity of 
the impurity. 

In summary, we have shown that the degenerate Anderson model is diagonalized exactly 
by the Bethe ansatz beyond LDA. If the impurity /-level energy lies inside the gap, the spec- 
trum of the charge excitations of the system exhibits a novel rich class of gap states, which 
have obviously been eliminated from previous studies due to LDA. It should be emphasized 
in conclusion that one of the most exciting possible applications of the approach developed 
could be an exact treatment of the superconductivity problem in the presence of magnetic 
impurities 0, where carrier dispersion cannot be linearized around the Fermi level of a host 
metal. 
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work. 
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